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Abstract 

We report on a study of exclusive radiative decays of the T(IS') resonance into the final states 
77r°7r°, ygg and yir°g, using 1.13 fb _1 of e + e~ annihilation data collected at y/s = 9.46 GeV 
with the CLEO III detector operating at the Cornell Electron Storage Ring. In the channel 
77r°7r°, we measure the branching ratio for the decay mode T(15) —> 7 / 2 ( 1270 ) to be (10.5 ± 
1.6 (stat) g (syst)) x 10 -5 . We place upper limits on the product branching ratios for the 
isoscalar resonances /o(1500) and /o(1710) for the 7r°7r° and gg decay channels. We also set an 
upper limit on the T(IS') radiative decay into 7r°?y. 

PACS numbers: 13.20.Gd, 12.39.Mk 
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Radiative decays of quarkonia, where one of the three gluons arising from the quark- 
antiquark annihilation is replaced by a photon leaving two gluons to form bound states, 
are thought to be a glue-rich environment that may lead to the production of glueballs and 
gluonic-mesonic states rather than ordinary mesons [l, 2j|. Lattice gauge theory calculations 
P, 0] predict that the lightest glueball should have J PC = 0 ++ and that its mass should 
be in the range of 1.45 to 1.75 GeV/c 2 , with decay into two pseudo-scalars ( J FC = 0” + ) 
expected to dominate. Unfortunately, the identification of a scalar glueball among the many 
established scalar resonances is difficult, as they have the same quantum numbers and similar 
decay modes and may mix. The triplet of /o states are likely candidates for the superposition 
of quark states and a scalar glueball state. Many of the lattice QCD models predict the 
decay ratios (e.g., 7777 / 7 T 7 t, rjrj/KK) for a glueball and for scalar resonances and this is a 
possible tool to distinguish among them. 

Most of the information on radiative decays of quarkonia has centered on J/-0 decays 
BBSS ilQj . leading to a list of two-body decay branching ratios. The establishment of a 
corresponding list for Y(IS') decays is desirable and would not only deepen our understanding 
of cc and bb quarkonia, but could also contribute to the identification of a scalar glueball 
state or shed new light on its mixing with ordinary nearby meson states. 

Recently, radiative decays into two charged particles have been studied by the CLEO 111 
collaboration [lT[ . The analysis included a measurement of the decay rate into /2(1270), a 
confirmation of its spin, and a measurement of its helicity distribution. In this analysis, we 
use the same CLEO 111 T(IS') data sample to perform a complementary study of all-neutral 
decays. Although these final states are subject to poorer resolution and efficiency than those 
with charged particles, they have the advantage of having no background from QED final 
states such as yp. Furthermore, they allow the search for states decaying into 7777 and 7i°r]. 
Resonant production in the latter mode would be a signature of unexpected physics. 

The analysis presented here uses data collected by the CLEO 111 detector configuration 
HQ at the Cornell Electron Storage Ring (CESR). The vital component for this analysis 


is the CsI(Tl) calorimeter, which has a resolution of 1.5%(2.2%) for 1 GeV(5 GeV) photons, 
typical of the photons studied here. We search for radiative Y(IS') decays in the modes 
T(IS') —■> 77 r° 7 r°, 77777 and 77 ^ 77 . The T(IS') data (E cm = 9.46 GeV) sample consists of an 
integrated luminosity of 1.13 fb -1 , corresponding to ( 21.2 ± 0.2(syst)) x 10 6 Y(IS') decays 


Candidate events for the individual final states ( 77 r° 7 r°, 77777 and n fK°r]) are selected in a 
similar fashion, using the following basic selection criteria. An event must have no charged 
tracks and exactly one electromagnetic shower in the barrel (| cosd| < 0.75, where 6 repre¬ 
sents the polar angle) or the endcap region (0.82 < | cosd| < 0.93) of the calorimeter with 
an energy exceeding 4 GeV, together with at least four other photons in the event. All 
combinations of two photons (excluding the photon that has E > 4 GeV) in the event are 
then combined to form 7 r° and 77 candidates. To be selected, an event must have two pairs 
of photons satisfying the requirement 


\M 2 K/ffi) + P%(A/ri2) < 5, 
with Pi and P 2 being the pulls, defined as: 

P( 7 r°/ 77 ) = m 77 — m(ir°/ri) /cr 77 , 


where m 77 is the 77 invariant mass, m(7i°/r 7 ) is the known 7 r° or 77 mass, 15] and <r 77 is the 
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77 mass resolution, with typical values of 5 - 7 MeV/c 2 . The 7r° and r) candidates are then 
kinematically constrained to their masses, m(7T°) and 777 ,( 77 )- 

To study the event-selection criteria and measure their efficiencies, we use a Monte Carlo 
simulation consisting of an event generator jl6j and a GEANT-based [ItJ detector-response 
simulation. For each final state, Y(IS') —> jX, events are generated with X = f 2 ( 1270), 
/o(1500) and / 0 (1710), using a Breit-Wigner line-shape and the PDG mass and width [15]. 
We do not search for the / 0 (1370) as it overlaps completely in mass with the / 2 (1270) due 
to its large instrinsic width. 

A 4-momentum cut and an asymmetry cut are then used to further select candidate 
events. For the 77 r° 7 r° final-state selection, the allowed region for the 4-momentum is 
bounded by the following three conditions: \p\ = —0.30 — 1.20 A E, \p\ = 0.25 — 0.80AF, and 
\p\ = 1.10 + 0.50 A E, where A E is the difference between the reconstructed event energy and 
the center-of-mass energy (E cm ) in GeV and \p\ is the magnitude of the reconstructed total 
event momentum in GeV/c. These cuts include the A E — \p\ area where the 4-momentum 
is conserved for the entire event and increase the efficiency by, in addition, including the 
region where the single, recoiling photon is reconstructed with too low an energy. For the 
latter, the 4-momentum is not conserved for the entire event but only for the intermediate 
resonance X in the decay chain T(IS') —> 7 A 7 " —> 77 r° 7 r°. These cuts are illustrated in Fig¬ 
ure [0 We define a 4-momentum allowed region for the 777 final state selection in a similar 
manner. 

A source of background originates from combining a wrong pair of photons to form a 7 r° 
or 7 candidate. Real 7 r° and 77 mesons decay isotropically and their angular distributions are 
flat. However, the 7T° and 7 candidates which originate from a wrong photon combination do 
not have a flat distribution in this variable and can largely be removed by a cut which uses 
the polar A 9 and azimuthal Atp angle differences between the two photons from a decay 
candidate. For the 77r°7r° final state the asymmetry requirement is yj A9 2 + Aip' 2 < 40°, 
while for the 777 final state it is ^A9‘ 2 + Aip' 2 < 60°. 
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FIG. 1: Illustration of the chosen 4-momentum distribution cuts on a larger (a) and smaller (b) 
scale, using Monte Carlo events for the decay channel Y(IS') —> 7 / 2 ( 1270 ) —> 77r°7r°. The slanted 
line in the lower right part of (b) divides the selected events roughly into two categories: the 
lower area where the 4-momentum is conserved for the entire event; the upper area where the 4- 
momentum is conserved for the intermediate resonance but not for the entire event. The remaining 
events outside the selected area in (a) have an energy loss from more than one photon and, hence, 
are excluded from the selection. 

Comparison of the invariant 7r°7r° and rjr) mass spectra from the Monte Carlo simulation 
reveals significant differences in the mass and width values from the ones used at the gener- 
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ator level. These differences are parametrized in the form of Gaussian resolution functions 
off-set from zero. The mass shift is an artifact of the shower reconstruction in the calorimeter 
of such fast 7 r° and 77 mesons, when the showers tend to overlap. For example, the resolution 
function for the decay /2(1270) —» 7 r + 7 r~ is a Gaussian function with a = 27 MeV/c 2 and 
an offset of —20 MeV/c 2 . 

We determine the selection efficiency for each of the resonances individually. The event 
selection efficiencies are summarized in Table 1. The uncertainties shown are statistical only. 


Resonance 

Reconstruction Efficiency in % 


77r°7r° 

77/77 

/ 2 (1270) 

16.4 ± 0.2 

10.2 ± 0.2 

/o(1500) 

20.4 ± 0.3 

9.1 ± 0.2 

/o(1710) 

20.6 ± 0.3 

8.6 ± 0.2 


TABLE I: Reconstruction efficiencies for various intermediate resonances in the 77 r° 7 r° and 77777 
final states. 


The major background contribution in our signal region originates from non-resonant 
processes. CLEO’s sample of data collected in the continuum below the T(IS') (192pb -1 ) 
is too small to perform a continuum subtraction. Hence, we parametrize the background 
using a threshold function of the form 

F(x) — N ■ (x — T) ■ e Cl{x ~ T) + C2(x " T)2 , 


where x is the 7 r° 7 r° invariant mass, N is a scale factor, T is the mass threshold, and < 7 , 
C 2 are free parameters. This functional form is a good fit to the spectrum obtained from a 
large Monte Carlo data sample of continuum events, and also to continuum events taken at 
energies near the Y(4S'). 

Figure 2 shows the final 7 r° 7 r° and 77/7 invariant mass spectra. The 7 r° 7 r° invariant mass 
distribution is dominated by the isoscalar resonance /2(1270). The 7777 invariant mass distri¬ 
bution, Figure 2(b), has only two events, which is too few to show any resonant structure. 

The Monte Carlo signal events for the processes T(15 l ) —> 7 X —» 77 r° 7 r °/7777 are produced 
with a decay angle distribution which is characteristic of the spin of the final state (he., J = 0 
for / 0 (1500) and J = 2 for / 2 ( 1270)). However, the generation does not take into account the 
correct helicity distribution for the / 2 ( 1270) since this distribution depends on the specific 
decay channel and can only be determined from the data itself. The method to obtain 
the correct helicity-angle distributions is described in detail in [ll[ and results in a helicity 
correction factor which takes into account decay-dependent efficiency corrections and the 
helicity substructure for the final state resonance. For this analysis, we use the helicity 


substructure, which is independent of the charge of the pions, determined in 11 ], as this 
is more precise than the one we can determine using the decay into 7 r° 7 r°. We obtain a 
correction factor for the /2(1270) of 0.66 ± 0.04, where the uncertainty is statistical only. 
This factor multiplies the efficiency stated in Table [7] 

To determine the branching ratio for Y(IS') —» 7/ 2 (1270), we fit the invariant 7 r° 7 r° mass 
distribution with a spin-2 Breit-Wigner line-shape of fixed mass and width, convolved with 
the resolution function derived from Monte Carlo studies as previously described, together 
with the threshold function. Integrating the Breit-Wigner line-shape fit from 0.28 to 3.0 
GeV/c 2 gives 67.9 ± 10.2 events for the / 2 ( 1270). 
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FIG. 2: The (a) 7 r° 7 r° invariant mass distribution, and (b) 7777 invariant mass distribution, from the 
Y(IS') data sample. The line shows the fit described in the text. 


With the results from this line-shape fit, the efficiency from Table Q] and the helicity- 
correction factor, we determine the product branching ratio for the / 2 ( 1270) to be: 

B(T(IS) -> 7 / 2 ( 1270 )) • B(f 2 ( 1270) ->• ttV) = (3.0 ± 0.5) x 10" 5 , 

where the error is statistical only. 

We determine a systematic uncertainty on this branching ratio of ±17%. The largest 
contribution to this error originates from uncertainties in the line-shape fit and the threshold 
function used for the background parametrization. Other contributions include systematic 
uncertainties in the 7T° reconstruction and in the 4-momentum cut. Taking into account the 
isoscalar nature of the / 2 (1270) and the branching ratio of £>(/ 2 ( 1270) 


7T7T ) 


= 0 . 847/°;°?2 


15], we determine an overall Y(15) radiative decay branching ratio to / 2 (1270) of: 


B(T(1S) -> 7/ 2 (1270)) = (10.5 ± 1.6 (stat) t\i (syst)) x 10 ~ 5 . 

To set upper limits on the branching ratios for other likely resonances in the 77 r° 7 r° 
final state, we include an additional spin-dependent Breit-Wigner line-shape in the / 2 ( 1270) 
branching ratio fit, with a line-shape determined from our Monte Carlo studies. We fix the 
area of the additional Breit-Wigner and then repeat the fit using different values for the 
number of events. We then plot the number of events versus their likelihood from the fit, 
numerically integrate the area under the curve and determine the number of events where 
90% of the physically allowed area is covered. This number represents the upper limit at 
the 90% confidence level (C.L.), which we find to be 6.9 events for the /o(1500) and to be 
6.6 events for the / 0 (1710). Using the branching ratio £>(/ 0 (1500) —> mr) = 0.349 ± 0.023 
15). and incorporating the systematic uncertainties (~6%) in the efficiencies by smearing 


the probability density function, we determine the 90% C.L. upper limit branching ratio for 
the / 0 (1500) to be 

B(T(IS) -> 7/o(1500)) < 1.5 x 10” 5 , 
and the product branching ratio for the / 0 (1710) to be 


B(?(1S) - 7/o(1710)) • £>(/o(1710) - ttV) < 1.4 x 10 


1—6 


As we see no evidence of any resonant structure in the 7777 invariant mass distribution we 
measure upper limit branching ratios for the /o(1500) and /o(1710). For this determination 
we use the simple method of event counting. The final invariant mass plot has negligible 
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background and, hence, we assume both events are from the Y(IS') —> 7 ?/?/ final state. 
Therefore, the number of events follows a Poisson distribution. For the /o(1500) we find 
1 event in the mass interval of 1 full-width around its mass and 0 events for the /o(1710), 
which translates into 90% C.L. upper limits of 3.9 and 2.3 events, respectively. 

The systematic uncertainty for the / 0 (1500) and the / 0 (1710) is ~ 30%. The largest 
contributions to these uncertainties originate from the 4-momentum requirement, and the 7 
asymmetry cut, which are on the order of 20%. Combining the statistical and systematical 
uncertainties, we determine the 90% C.L. upper limit on the product branching ratio for the 
/o(1500) to be: 


B{?{1S) - 7/o(1500)) ■ £>(/o(1500) - 777 ) < 3.0 x 10~ 6 , 
and for the / 0 (1710) to be: 

B(T(15) - 7/o(1710)) ■ £>(/o(1710) - 777 ) < 1.8 x 10- 6 . 

In the decay Y(IS') —> 7 A", if we assume that the 7 is produced directly and is not 
the product of an intermediate virtual particle, the resonance X must be an iso-scalar. In 
this case, if X is conventional meson state, it can only decay into a pair of pseudo-scalars 
( J p = 0") each with 7 = 0 (e.g., 77 ), or I = 1, e.g., mr. Observation of a resonance in 
7 t°7] could therefore be an indication that the photon in this case is the result of enhanced 
production via an intermediate hadron, or alternatively the result of an unexpectedly large 
1=0 component of the 7 r °rj final state. 

Following the same analysis chain as detailed above and using the same p, A E region as 
for the 7 r° 7 r° case, we find no events in our signal region for this decay. Hence, we determine 
an upper limit for the branching ratio Y( 1 S) —> 77 r 0 ?/. 

We use the same method as for the upper limit determination in the Y(IS') —> 7777 final 
state. To measure the reconstruction efficiency for any exotic-state mass, we generate Monte 
Carlo events of the type Y(IS') —> 77^7 with a flat 7i 0 r] invariant mass distribution between 
0.7 and 3 GeV/c 2 , and use the lowest efficiency found in the entire mass distribution of (4.8 
± 0.5)%; The efficiency is relatively flat over the mass interval of interest. 

Having no events in the data over the mass range of 0.7 to 3.0 GeV/c 2 corresponds to a 
90% C.L. upper limit of 2.3 events. Combining this with a systematic error of tu%, due to 
the same sources of uncertainty as with the previous 2 analyses, we determine the 90% C.L. 
upper limit for the branching ratio to be: 

B{T{IS) ^ 7 A) < 2.4 x 10 ~ 6 . 

In summary, we have analyzed 1.13 fb -1 of data from the CLEO III detector at the Y(IS') 
for resonances in the radiative decay channels Y(IS') —> 77 r° 7 r°, 7777 and 'yir 0 r]. 

In the decay channel 77r°7r°, we measure a branching ratio value for the isoscalar resonance 
/ 2 ( 1270) of B(T(1S) —> 7 / 2 ( 1270 )) = (10.5 ±1.6 /};§) x 10 -5 . This is in excellent agreement 
with the same branching ratio obtained from the charged final state 77r _ 7r + , using the same 
CLEO III data set: B(T(1S) —> 7/ 2 (1270)) = (10.2 ±0.8 ±0.7) x 10~ 5 jll[. It also agrees 
within the uncertainties with the earlier CLEO II result of (8.1 ± 2.3 ± 2.7) x 10~ 5 , based 
on the decay channel 77r _ 7r + [18j; this earlier measurement had no correction for the hclicity 
distribution, and the large systematic uncertainty reflected this fact. 

In addition, we determine 90% C.L. upper limits for the isoscalar resonances /o(1500) 
and /o(1710) decaying into tht, as well as a 90% C.L. upper limit for the decay Y(IS') —»■ 
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7 /o (1500). Based on the scalar-glueball mixing matrix from [5], QCD factorization model 
calculations in [2] predict branching ratios for the /o(1500) and /o(1710) to be £>(Y(1S') — > 
7/o(1500)) « 42 —84x 10 -5 and i3(T(lS') -> 7/ 0 (1710))-£(/ 0 (1710) -»■ 7r°7r°) « 6-12xl(T 6 . 
Our measurements of B(T(1S) —> 7/ 0 (1500)) < 1.5 x 10 -5 and £>(T(1S') — > 7/ 0 (1710)) • 
S(/o(1710) — > 7r°7T 0 ) < 1.4 x 10 -6 are much smaller than these predictions. 

In the 7777 decay channel, no resonant structures are observed. Therefore, we determine 
a 90% C.L. upper limit on the branching ratios for the isoscalar resonances /o(1500) and 
/o(1710) decaying into 77 as £>(Y(1S') — ► 7/ 0 (1500)) ■ £>(/ 0 ( 1500) —> 77 / 7 ) < 3.0 x 1CT 6 and 
B{?(1S) -> 7/o(1710)) • 8(/ 0 (1710) -»• rjrj) < 1.8 x 10 - 6 . 

The search for states in the 77 r °?7 decay channel does not show any evidence of a signal. 
We determine a 90% C.L. upper limit on the branching ratio for the decay Y(IS') —> 77 r°r/ for 
any intermediate state with a mass between 0.7 and 3.0 GeV/c 2 to be /5(T(1S') —> 77 r°? 7 ) < 
2.4 x lO" 6 . 
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